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Possible Revelation of Seesaw Mass Pattern 
in Solar and Atmospheric Neutrino Data 
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Abstract 

Assuming the solar and atmospheric neutrino deficits to be due to neutrino oscilla- 
tions, it is shown that the 3x3 mass matrix spanning v e ,u^, and v T may have already 
revealed a seesaw mass pattern. Also, this matrix is the natural reduction of a simple 
5x5 seesaw mass matrix with one large scale, the 4x4 reduction of which predicts 
that a fourth neutrino would mix with v e and in such a way that i/„ — > v e oscilla- 
tions may occur just within the detection capability of the LSND (Liquid Scintillator 
Neutrino Detector) experiment. 



Whether or not neutrinos have mass is clearly a fundamental issue in particle physics and 
astrophysics. There is now a good deal of indirect evidence from measurements of the solar 
v e flux[|I], [1 and the atmospheric v^jVf, ratio[^|, |5| that neutrinos oscillate from one to 
another. This means that v e and are not mass eigenstates; each is rather an admixture 
of two or more mass eigenstates. Under the simplifying assumption that only two neutrinos 
are involved in each case, probable central values for the mass-squared difference and the 
mixing angle are 

Am 2 ~ 4 x 10" 6 eV 2 , sin 2 29 ~ 0.01, (1) 
from analyzing]?]] solar data, assuming the small-angle, nonadiabatic solution, and 

Am 2 ~ 10" 2 eV 2 , sin 2 29 ~ 0.5, (2) 

from analyzing || atmospheric data. 

Given only the above information, it is clearly not possible to reconstruct unambiguously 
the underlying 3x3 neutrino mass matrix. However, if a certain empirical relationship can 
be found among the four parameters listed above, an important insight may be gained as to 
the form of this mass matrix. Such an example is already very well-known in the case of the 
quark mass matrix. It was pointed out many years ago|| that the empirical relationship 

sin 2 9 C ~ m d /m s (3) 

may be obtained with a 2 x 2 mass matrix of the form 



M 



a 
a h 



(4) 



This simple observation has generated over the years an enormous literature on quark mass 
matrices. It is an especially active field of research in the past two or three years. In the 
jargon of neutrino physics, the form of M. in Eq. (4) is called seesaw [[K] . 
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We propose here that Eq. (4) is also applicable to the 2x2 matrix spanning v T and the 
linear combination cos 9 + v e sin 9, where 9 is the angle measured in atmospheric neutrino 
oscillations. The orthogonal linear combination v e cos 9 — sin 9 is assumed to be massless. 
The 3x3 neutrino mass matrix spanning v e ,v^ v T is then given by 



M, 



b sin 2 9 b sin 9 cos 9 a sin 9 
b sin 9 cos 9 b cos 2 9 a cos 9 
a sin 9 a cos 9 



(5) 



The eigenvalues of M. u are 0, —a 2 /b, and b for a «b. Let the corresponding mass eigen- 
states be U\, u 2 , and v%. Then the mixing matrix U is easily obtained: 



V »r ) 



I cos 9 — (a/b) sin 9 sin 9 \ 

— sin# —(a/b) cos 9 cos9 
V 1 a/b 











) 


V ^3 / 



(6) 



Following the results of the two-neutrino analyses given in Eqs. (1) and (2), let b = 0.1 eV 
and a = 0.014 eV, then a 2 /b ~ 2 x 10~ 3 eV, hence Am 2 2 ~ 4 x 10~ 6 eV 2 and Am 2 3 = 
10 -2 eV 2 . Now choose sin 2 29 = 0.5 as given in Eq. (2), then the effective sin 2 29 12 = 
4\U e 2\ 2 \U e i\ 2 /[\U e i\ 2 + \U e2 \ 2 } 2 for solar neutrino oscillations is predicted to be 



sin 2 2#i2 ~ (a/b) 2 sin 2 29/ cos 4 9 ~ 0.014, 



(7) 



which roughly matches the value given in Eq. (1). This marks the first time that a possible 
seesaw mass pattern in neutrino physics has been identified. 

Since Amf 2 << AmJ 3 and a/b << 1, the atmospheric neutrino deficit |5|, [J is explained 
by a simple two-neutrino oscillation between and v e . However the solar neutrino deficit 
comes from the oscillations of all three neutrinos. If matter effects fll[ could be neglected, 
the disappearance probability of solar neutrinos would be given by 



P{v P 



2\U^\ 2 \U el \ 2 + 2\XJ eZ \ 2 \U e2 \ 2 + 2\XJ e2 \ 2 \XJ el \ 



cos ■ 



t/S.m\ 2 
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where U a i are the elements of the mixing matrix in Eq. (6). We see that the effects of the 
heaviest mass parameter average out, but a dependence on \U e 3\ 2 remains. Only in the limit 
of vanishing \U e s\ 2 = sin 2 9 is the two- neutrino approximation valid for solar neutrinos. This 
general observation is also true when matter effects are included |T2|. Since a large sin 2 # is 
required by the atmospheric data, the naive comparison of Eq. (7) with Eq. (1) will have to 
be modified. 

A three-neutrino analysis of all the solar neutrino data|l], || f|, f§ is shown in Figs, 
(la), (lb), and (lc). The method is the same as that of an earlier analysis [|T3p where a 
continuous range of \U e ^\ 2 values is considered. Three values of Am 2 3 and sin 2 29 which 
fit the atmospheric neutrino data have been chosen, and their predictions for Am 2 2 and 
sin 2 29 e2 = 4|f/ e2 | 2 (l — |£/ e 2| 2 ) (dashed line) have been compared to the solar data constraints 
(solid contours). For Fig. (la), a small value of sin 2 29 is chosen. There the allowed solar 
neutrino solutions are close to those of the two-neutrino analysis, and the predictions of 
Eq. (5) match the small-angle, nonadiabatic solution. For Figs, (lb) and (lc), two larger 
values of sin 2 29 are chosen and the allowed solar neutrino solutions are very different from 
those found in the two-neutrino analysis. Here the predictions of Eq. (5) match an adiabatic 
solution which is allowed only in the three- neutrino analysis. As shown in Fig. (2), Eqs. 
(5) and (6) are satisfied in different ways across the parameter region indicated by the 
atmospheric data. Although we do not include the region allowed by the new Kamiokande 



multi-GeV data|| which is in conflict with Frejus [p!4 | , solutions exist there as well. 



Consider now the possible origin of the 3x3 neutrino mass matrix hA v of Eq. (5). The 



most natural way is to find it as the seesaw |T0| reduction of a larger matrix which includes 
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a higher mass scale. A very simple solution is the following 5x5 mass matrix 












mi 














m 2 

















m 3 


mi 


m 2 








m 4 








m 3 


m 4 


m 5 



(9) 



spanning u e , v T , u s , and N, where the last two are singlet Majorana neutrinos. In the 
limit of very large m 5 , the seesaw reduction of M. 5 is 



Ma = 



mi 

m 2 

—ml/m 5 ~m 3 m 4: /m 5 



(10) 



mi m 2 — 77137774/7715 — 777I/7775 
Assume now that 7774/7775 is the dominant term in M.^ then a second seesaw reduction gives 
the following 3x3 matrix: 



m\m^jm\ m 1 m 2 m 5 /m'l —mim^/m^ 
mim 2 m^jm\ m\m^jm\ —m 2 m^/m4 
—777! 7773/7774 —77727773/7774 



;n) 



which matches exactly M. v of Eq. (5) with the identification 



b = 



(m\ + 7772)7775 
ml 



tan 9 = 



mi 



7773 



a = 



m\ 



(12) 

7772 77*4 

With the further requirement that 77737774 << m^\jm\ + m\ so that a « b, we then obtain 
Eq. (6). Since Mz(= M u ) reduces further to M of Eq. (4), one might even speculate that 
the present neutrino data are indicative of three successive seesaw reductions. 

If we take .M5 (and thus seriously, then we have one more very interesting prediction. 
There should be a fourth neutrino (1/5) with mass given by —m\jm^ and mixing to v e and 
given by — 777 1 m 5 /m\ and —m 2 m^jm\ respectively. Since 



m\\ /mim 5 \ I m 2 m^\ f{m\ + m\)m^ 



7775 J \ ml J 



m 1 



7774 




777i 



ml + 7772 




777 2 



mf + 7772 , 



(13) 



we predict that there should be additional oscillations between v e and through this fourth 
neutrino mass eigenstate with restricted values of Am 2 and sin 2 29 such that their product 
is the same as that for atmospheric oscillations. Using the numerical results obtained earlier 
in this paper, we would then predict 

3 x KT 3 eV 2 < Am 2 sin 2 29 < KT 2 eV 2 . (14) 

Part of this region is excluded by the E776 neutrino experiment at Brookhaven National 
Laboratory ||16|| , but another part (Am 2 ~ 1 eV 2 , sin 2 29 ~ 5 x 10~ 3 ) lies within the detection 
capability of the LSND (Liquid Scintillator Neutrino Detector) experiment at Los Alamos 
National Laboratory. We also note that the upper bound in the above comes from reactor 
experiments ||15|| . Relaxing it would allow a larger Am 2 (~ 6 eV 2 ) which is preferred by the 
preliminary LSND results, but in conflict with the published E776 data. 

To obtain Ai$, the following softly broken discrete Z 3 symmetry may be considered. Let 
z/ e , z/^, z/ T , us, and N transform as uo, lu, 1, uj 2 , and 1 respectively, where a; 3 = 1. Then mi, 1712, 
and ?»3 come from the vacuum expectation value of the standard Higgs doublet, and m 5 is 
allowed by Z 3 . However, m 4 breaks Z 3 softly and so does the diagonal mass term for ug, 
but as long as the latter is much smaller than m 2 /ms, the reduction to Ai u proceeds as 
before. With four light neutrinos, the nucleosynthesis bound|T7| of N u < 3.3 is an impor- 



tant constraint. Although v$ is a singlet neutrino, it mixes with other neutrinos and may 



contribute significantly to N v through oscillations |18|. In the context of .M4, this constraint 



implies that y mf + m 2 . has to be very much smaller than m 2 /ms, hence the sin 2 29 param- 
eter in Eq. (14) is too small to be in the range of the LSND detection capability. In other 
words, the simplest extension of our basic proposal to four neutrinos cannot accommodate 
simultaneously the LSND results and nucleosynthesis. The latter constraint can be avoided 
if we switch v T and z/g in but then the underlying theory becomes much more involved. 
For example, the v e v T and f^v T mass terms must now come from either a Higgs triplet or a 
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radiative mechanism as recently proposed]!^]. Details will be given elsewhere. 

In conclusion, we have shown in the above how present data on solar and atmospheric 
neutrinos may have revealed a seesaw mass pattern for the three known neutrinos. This 
seesaw scenario requires the solar neutrino deficit to involve the oscillations of all three 
neutrinos. Hence a much more extensive numerical analysis is needed |l3j beyond that of the 
usual two-neutrino assumptionQ. It is interesting to note that with only two neutrinos, the 
adiabatic solution is not allowed by the combined solar data, but with three neutrinos, it 
becomes allowed if the mixing of the third neutrino with v e is large enough. 

We have also shown that a simple 5x5 mass matrix with one large scale reduces naturally 
to the desired 3x3 M. v of Eq. (5). Furthermore, the intermediate 4x4 reduction predicts 
a fourth neutrino which mixes with v e and in such a way that — > v e oscillations may 
occur just within the detection capability of the LSND experiment. 
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Figure Captions 



Figs. (1) Plots of Affl| 2 versus sin 2 29 e2 = 4|6 r e 2| 2 (l — |^ e 2| 2 ) for various values of sin 2 26 l . 
The solid contours surround the parameter region allowed by the solar neutrino data 
EL H> H' M & t 90% confidence level. The dashed line gives the prediction of Eq. (5) for 
a specific value of Am 2 3 which explains the atmospheric neutrino deficit. 

(a) sin 2 29 = 0.35 and Am 2 3 = 3.0 x 10~ 2 eV 2 . 

(b) sin 2 29 = 0.50 and Am 2 3 = 1.0 x 10~ 2 eV 2 . 

(c) sin 2 29 = 0.75 and Am 2 3 = 4.0 x 10~ 3 eV 2 . 

Fig. (2) Plot of Am 2 3 versus sin 2 29. The dashed contours surround the parameter region 
allowed by sub-GeV atmospheric neutrino measurements |6|] but excluded by the Frejus 
data[TJJ] and below constraints from reactor measurements flTBf . The shaded regions are 



where Eq. (5) can satisfy the solar neutrino data; the left shaded region corresponds to 
the small-angle nonadiabatic solution and the right shaded region corresponds to the 
adiabatic solution. In between these two regions is a small area where it is not possible 
to satisfy the solar neutrino data. 
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